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ABSTRACT 


Suspended load is continuously deposited on the banks of 
alluvial rivers with active channels. An increase in suspended 
PCodgwiLhoug chance in hydraulic conditions is, thought to 
induce narrower channels. Stream-bed instability of laterally 
Pranew beds initiaved by ithe effect. of suspended..load stransport 
in the transverse direction, deepening the centre and building 
the sides, can explain the geomorphological phenomenon of bank 
formation. 

ALlinear anaiysis of lateral instability of an .erodible 
stream-bed has been carried out. Forcertain conditions a 
growth in the amplitude of bank-like perturbations confirms the 
existence of stream—béed instability in the lateral direction 
resulting in the formation of banks. The question as to why 
Spontaneous bank formation is not observed in laboratory flumes 
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CHAPTER: 1 


INTRODUCTION 


1.1 GENERAL STATEMENT 

Knowledge of the instability of the fluid-bed interface 
has been extended in recent years, allowing for investigation 
and explanation of various aspects of the geomorphological 
behavior of natural streams. Interaction of flowing water and 
an erodible bed results in various bedforms. Many of the same 
methods used to analyze such bedforms as dunes can be used to 
explain the formation of banks, which are actually large scale 
bateralsbed-configurations. In fact, equilibrium banks can be 
formed in overly wide channels as a result of an erosion- 
deposition process, caused by a systematic perturbation of the 


ross: lateral transport of sediment. 


ie OBJECTIVE 
The objective of this analysis is to investigate the 
process of bank formation. The formulation of A ere nace 
model for this purpose. involves five essential steps: 
1. The equattons: of sediment conservation describing the 
erosion-deposition process, subject to the limitations 
imposed by two boundary conditions, i.e. those at the 


water-surface and the bed, are formulated. 
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The equations are linearized for slight perturbations 
about steady ant oun flow. 

The equations are put in non-dimensional form which 
Drovides'-an appropriate format for. the final results. 
sinusoidal perturbations characteristic of incipient 
bank formation are introduced and the associated dis- 
persion equation is analyzed. 

Stable and unstable regimes of initially laterally flat 
ehannels are established. The latter corresponds to 
incipient bank formation. The associated conclusions 


constitute the final stage of the analysis. 
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CHAPTER 2 


REVIEW OF LITERATURE 


2nd DEE IN ETIONS 

Some of the terms frequently used in this report are 
defined below. 
ins tabiLity 

If the equilibrium state of a system, when perturbed 
SHightly, exhibits a growth of initial perturbation, that equi- 
librium state of the system is said to be unstable. | 
Turbulent Diffusion 

Turbulent diffusion is the process by which random motions 
satisfying the momentum and mass balance of fluids spread mass, 
momentum, etc. from regions ofr nigh concentration to regions of 
low concentration. 
Bed Load and Suspended Load 

Bed load is the relatively coarse part of the sediment load 
that moves along the stream bed. Suspended load is the rela- 
tively fine part of the sediment load that is distributed 
throughout the wertical section of a stream. 


Erosion and Sedimentation 


Erosion is the removal of soil particles from their environ- 
ment. Erosion in streams refers to the removal of sediments by 


water. Sedimentation refers to the deposition of eroded soil 
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Erosion and sedimentation, as referred to in this analysis, 
ere sconrined: tO noncohesive material only. The reason for this 
is that the behavior of cohesive material is different from 
that of noncohesive materials: its complicated physiochemical 
properties would make the study intractable. 

Dunes and Antidunes 

The term dunes is commonly used for bed features whose 
VertaGpion in elevyauion is more “or less out of phase with that 
of surface wave above, while antidunes are more or less in phase 
with surface elevation. Dunes and antidunes are shown in 
Pugure 2.1: 

Ripples 

individual dune patterns are referred to asjripples which 

are understood to be small bed forms, whereas dunes are large 


ones. 


2.2 EROSION AND SEDIMENTATION IN NATURAL STREAMS 

2.e-l1 Conmecept 

There are two important categories of channel erosion: 

1. Stream-bed erosion 

2. Stream-bank erosion 

While bed erosion is mainly in the direction of flow, bank 
erosion is a lateral process, converted to a downstream pro- 
cess only as the eroded bank material is carried downstream 
or causes sedimentation on the stream-bed. Shear stress governs 
the erosion process in the longitudinal direction of flow. 


However, on a channel bank another force plays an equally 
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important role - the gravity force that causes particles to move 
down the sloping w#ide of the channel. Lateral transport of 
suspended matter gives rise to sedimentation. It is the inter- 
action between sedimentation and erosion that is assumed to 


form banks. 


2.2.2 Measuring Erosion and Sedimentation 
For channels in noncohesive sediments a tool for qualita- 
tive prediction of longitudinal erosive channel conditions is 


provided by Lane's (1955) relationship: 


QS « a G.. qd. 
where, Q .= stream discharge 
S = longitudinal stream slope 
G. = bedvsediment discharge 
d., =a characteristic diameter of the bed material 


ini proportion is particularly useful when two of the four 
Variables can be assumed to remain constant. 

Quantitative estimates of gross channel erosion or sedi- 
mentation are obtained from time sequence comparisons of surveyed 
cross-sections (Vanoni, 1975). They can be indirectly estimated 
by relating sediment discharge to flow regime, a change in 
which can reflect the rate of erosion of sedimentation. 

The above measurement techniques may be applied to a river 
that exhibits»an unstable section due to excessive sedimentation 
Or erosion: but formarriverminatcis statistically stable, i.e. 


a river that is neither widening nor narrowing its channel on 
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an average basis, very little information is available to 


compute either the rate of “erosion or’ deposition of bank material. 


2,5 RIVER SEDIMENT; MODE OF TRANSPORT 

Zu Sek.) COnGen tT 

DW Senerat wana river, sediments raneineg in size from clay 
particles to boulders may be found. Depending upon the mode of 
transportation, the sediments are classified as either bed load 
Oreo uspended 1oad, “As far as sdze is concerned, there is no 
elear demarcation between suspended and bed load. A particular 
Fiver under 4 particular flow situation may transport specified 
Size ranges of suspended and bed load materials; but this range 
may not correspond to another river with a different flow condi- 


GLOn:. 


2.3027 oCdimens Load ‘Transport Theertes 

The mode*of transport of bed load ts *different from’ that of 
suspended load; the former being transported in a saltating, 
or jumping mode while in close proximity to stream-bed; and 
the latter being transported in suspension while being distrib- 
uted over the entire section of the stream. The precise 
mechanism *of'bed load “transport' ts “out-of the scope of this 
study. The theory of suspended load transport is described 


separately in the following chapter. 


2.4. THE THEORY OF SUSPENDED LOAD 
2.4.1 The Turbulence Models 
The suspension of sediment particles in a flow is attributed 


to turbulence. The turbulence model as formulated by Kennedy 
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(1895) states that the whole body of water is kept in circulation 


que to turbulence; this favours suspension. 


Some of the earlier studies of turbulent suspension reveal 


three important facts that are also supported by experiment. 


de 


As regards_to the vertical distribution of suspended 
matter, more suspended particles move near the bottom 
than near cme water surface. This is, illustrated in 
Figure 2.2, in which the suspended sand concentration 
distribution of four different sand diameters are 
plotted. The data are due to Kalinske and Pien (1943), 
col¥ected at the centre of a 2.25 foot wide daboratory 
flume. 

Thé lateral distribution of suspended matter’is such 
that the Suspended sediment concentration® is higher 
near) the centre than near the banks This a2 iidus— 
Crated inericure 203, Waren is produced fromecher data 
of Hubbell and Matejka (1959), collected at three 
different times during one year period near Dunning, 
Nebraska. 

A systematic variation of suspended sediment concentra= 
tion in the downstream direction: is probabiyenot present 


over fairly short, uniform reaches#with no tributaries. 


2.u%2 Diffusion —Dispersion Models 
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The suspension of sediment particles due to turbulence is 


treated in the model as a diffusion-dispersion process. In this 
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process, the particles are diffused in the medium while being 


transported. 


2.4.2.2 Equations of Susvended Sediment in Shear Flow 

In a laminar flow, diffusion means the spreading by random 
mMOctLON; “which o¢ecurstwhether the fluid) is moving or not. 

it*c isthe concentration’ ofa conservative substance’and 
re is the molecular kinematic diffusivity, the equation of 


eonservation- of thersubstance in a motionless fluid is given by, 


ac = Hey=e (2.19 


dt 
Bduavion 2.1 deseribes the diffusion process. If.the velocity 
vecvor of the laminar flow is OL 


ac 
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Equation’ 2.2 describes the dispersion process, which combines 
the action of convection and diffusion. ~in a turbulent Plow, 
molecular diffusion,is, far.weakerrthan the turbulent, diffusion 
induced by convective motion of the eddies. Thus in an averaged 
form of (2, 2)ethe scalar molecular, diffusivity can be replaced 
With wnef fective eddy diffusivity tensor. <Atethis, tensor “is 
assumed to be diagonal, the mean three-dimensional equation for 


the conservation of suspended material takes the form, 


ee ,y Mey May q cy (2.3) 
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where, c and (Uy, Uy» Uy) are now interpreted to be mean quantities. 
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Here bx? by» and eo are the turbulent kinematic eddy 
Guerusivities in x;y, and z directions respectively. Also here 
xX, Y, Z correspond to longitudinal, lateral, and vertical 


direction coordinates respectively. 


Deore  Nertical Eddy Difrfusivi ty 
rie vertical eddy diffusivity of flutd mass aT and the 
vertical eddy diffusivity of suspended mass E> are piven by 2 


general relationship, 


where, *m is) 2 CONStant ol propertionalityr Grat. (1971) cites 
fromthe literature *that ‘the value of m can be approximated as 
sige 

If mis the vertical turbulent mixing coefficient of 


momentum, the form due to Rouse (1937) can be written as, 


spa Z. ( 
m= 5 (L- 5) KDU, i225) 
where, K = Von Karman's constant 
DD t=) Depth of flow ins2 direction 


U, = Shear velocity at the bed 
Experimental data from Kalinske et al. (1943) suggest that 
the turbulent transfer coefficients of mass and momentum are 


approxi imavely equals a .e:. 
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Hence, we may approximate that E> = by so that ED can be given 
by, 


par Zz (2 eSu) 
g hott RDU 
ane distribution of ay is sec sO (De parabolic in the vertical 
Grrection, a6 reperteds by Kalinske et al... This tis il@ustrated 
in Figure 2.4. 

Fischer (1973) suggests that for practical purposes, the 


depth-averaged value of ae 
c=) 087 DU, C2.7) 


May.be wsed in plage of (2.6)... 


Engelund (1970) suggests the value, 
E = .0G DU, 2.8) 


2.4.2.4 Transverse Turbulent Eddy Diffusivity 

The nonisotropic mature of turbulence makes it difficult 
to establish an analytical relationship between vertical and 
transverse eddy diffusivities. Fischer (1973) provides a de= 
baaved discussion of empirical ames of the transmerse eddy 
diffusivity determined from experiments. These experiments are 
divided into three categories: (1) laboratory experiments 
with floating particles, (2) laboratory experiments with dis- 
solved tracers, and (3) field experiments with dissolved tracers. 

The laboratory floating particle experiments as mentioned 


above, yield an average surface dimensionless mixing coefficient 
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of by/ DU. = 0.20. An average of all laboratory experiments 
with dissolved tracers yields a ddvailatebedes transverse mixing 
coefficient Be) Be = 0.15. In curving channels, the transverse 
mixing coefficient can be greatly accelerated by secondary 
eurrents > ’'Pischer's findings indicate that -in terms of ‘DU,, 
the transverse mixing coefficients for actual channels range 
from Ey/ Us Pay 0 abatato.c2). 1, 

Fisher quotes Okoye (1970) in his correlation of oy 
against the aspect ratio B/D, where B is the width of the 
channel. For values of B/D greater than 10, the depth-averaged 


value of by/ Us is»seen to be between 0.1 to 0.22. 


2.4.2.5 Vertical Distribution of Suspended Sediments 
Assuming a steady-state distribution of suspended matter 
Ehave1s Unitormign the Tongeitudinal and transverse directions, 


equation 2.3 can be integrated to yield, 


: de 
QO = wae a b. ae (2.9) 


where, ie is thesparticle settling velocicy . 
mntecrating tequation 2.9.5) the solution is found to pe, 
- te (z-a) 
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Z 
c/Ca = e C25UO) 


Here ao tse) taken as’) a constant and Ca is ‘a’ reference ‘concentra— 
tiontiaies hHeient Sate trom the ped.Aelf a vertically constant 


value of oo isachosent, ‘a’ can be quoted to zero. 
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Equation 2.10 indicates that the concentration distribution 
is an exponential one, the concentration being lamegery.c Losier bo 
the bed than farther away from it. The steepness of the expo- 
nential curve depends upon the coefficient of diffusivity as 
and the settling velocity Vee Comparison with field data, 
mainly for wide rivers, leads Lane and Kalinske (1941) to the 
conclusion that, Equation 2.10, as approximate as it may be, is 


Ssuificiently accurate for practical purposes. 


Engelund (1970) has found that with the use of equation 2.8, 


the nominal bed concentration Ca for a = 0 can be represented 


as, 
Co = 0.0073 \Us 3 uy 
ive 
Ss 
Ux 
This equation is acecurate only,up to cone xy Hor. -vaLues not 
Ss 


a8 greater than 3, the volume concentration at the bed level 
Ss 


becomes ilarger,-than 0.32. -Thisavalue-is suggested.to be the 
maximum attainable volume concentration on the basis of experi- 
mental studies in flumes and field measurements of hyperconcen- 
tration (see section 2.5). 

Another approach to suspended load distribution is due to 
Rouse (1964), who solves equation 2.9 with the form for Eo 


reported in the equation 2.5 to find, 


Zt 
10 tecne ) 


Where. dh as (Me , is the Rouse Number. 
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pe reavure- er equation’ 2.13°is°that’z°= O° (bed level)? the 
concentration becomes infinite; and hence the equation does not 


apply right at the bed where z= 0. 


2.5 HYPERCONCENTRATION OF SUSPENDED SEDIMENT 

eno .k, Concent 

Beverage et_al. (1964) have defined Eecr eon entestions to 
be those that exceed 40 percent by weight or 15 percent by 
volume. 

The occurrence of extremely high, suspended sediment 
concentration is common to several alluvial streams. Beverage 
cv ol. Have compiled’ some examples or @hyperconeéentrated-reaches 


of’ such streams. 


Povwe Hyperconcentration and Turbulence 


The .etfect of high concentration on turbulence has been 
studied by Bagnold (1955), Nordin (1936) and Hino (1963). Many 
authors have discounted the role of turbulence as the major 
Factor in suspending high! concentrations, of sediment. High 
damping of turbulence due to increases in suspended sediment 
concentration has been suggested by Bagnold. On the other hand, 
Hrmo predicts that the decrease in turbulence Antensity «due co 
hyperconcentration of suspended material should be very small. 
Hino's conclusion results from an analysis of the energy 
equation of flow with suspended particles, taking into account 
the major factors liable to influence flow with suspension. 


Beverage et al. quote from the literature that even at 68 percent 
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concentration by weight (28 percent by volume), streams have 
been observed to possess all the characteristics of fully- 
developed turbulence, with numerous intermittent sand waves and 
high velocities. In any case, the amount of damping realized 
in hyperconcentrations appears to be strongly dependent on 


penercke usize: 


2.o058 9 Hy¥perconcentrationvand the Von. Karman Constant 

In open channel flow without sediment, the Von Karman 
Constant Ki i.e. the universal constant in the logarithmic 
velocity law, is given approximately by K = 0.4. With an in- 
erease .in suspended sediment concentration, K decreases. Graf 
(1971) quotes from the literature stating that a reduction in 
K-rseduertostherdamping effectyor eturbulence sichinsteim étidal . 
(1954) plot the nate gofeexpenditunesofrinictional (energy ines 
quired Go'support ithe suspendedrpartickes versus *K.; The graph 
indicates that-K decreases with increase in the rate of energy 
expenditure pewhieh tisediréecthyiproporvionak toetheraverage 
concentration of suspended sediment. With the suspension of 
neutrally buoyant particles in water, Hino's investigation 
indicates that K decreases, and the turbulent intensity increases, 


with an increase in volume concentration of particles. 


226 STABILITY OF ERODIBLE BEDS 

2.0)..5. COmeert 

The mathematical approach of Helmholtz (1888), quoted by 
Graf (1971) and Lamb (1945), indicates that a boundary between 


two fluids of different densities moving with different velocities 
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ts fsubject.to instability leading to,interfacial)waves. Graf 
quotes Baschin (1899) who suggested that the Helmholtz insta- 
bility could explain bedforms (see chapter 2.6.2) if the loose 
sand can be considered to act like a fluid. Kennedy (1963) 
indicates that the instability of the sediment-fluid interface 
depends upon the depth and velocity of flow and the properties 
of sediment and fluid. “The origin of meandering and braiding 
of alluvial rivers has been analysed using stability theory 
by Engelund and Hansen (1967), Callander (1969) and Parker 


(1976), to mention a few. 


2.6.2 Flow Regime and Bedforms 


In open channel flow the Froud number NF is often used as 
a-tlow eriterion. «For the purpose of classification of bed- 
forms, three flow regimes may be distinguished: (1) a lower 
regime with NF < 1; (2) an upper regime with NF > 1; and (3) a 
transition between the upper and lower regimes with NF = l. 

Idealized sketches of various bedforms are shown in 
‘Figure 21+. The definitions of- ripples, idunes, antidunes have 
been provided earlier in Section 2.1. Ripples and dunes usually 
occur in lower flow regime while antidunes are confined to upper 


flow regime. Washed-out dunes occur in the transition zone of 


the flow regime. 


2.6.3. Model Analysis 


22643. 1) Concent 


Stability analysis of the stream-bed interface can provide 


theoretical predictions of the geometry and behavior of bedforms. 
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Different attempts to explain bedform mechanics have resulted 
in nydraulic,.povential arid rotational model analyses, some 
involving the stability analysis of the fluid-bed interface and 
some Using the classical hydraulic approach. Reynolds (1975) 


provides a detailed description of each of the model studies. 


2.6.3.2 Hydraulic Models 
Graf (1971) quotes Exner's (1925) analytical treatment of 


an erodible stream-bed using a one dimensional flow approach. 
Exner’s model is classical and is the predecessor of the 
hydraulic model. jin this model} velocity variation with depth 
is reduced to a sectionally averaged mean velocity. 

ReLerpringe To rigure. 2.5. the) governing equation in Exneéer's 


mode lwis the .oquationsoh cote rnuicy.swhich may be written, as, 
(ii =) iB Ga) NE ay ROS CONS Cans 62513) 


where, h is the watersurface height 

mi is: the bed heient) tromvan arbivrary, datum 
U is the sectionally averaged velocity 

B(x) is the channel width that varies in the x direction 

of flow 
QO is thesflow rate 
In progressive stages, the model takes into account effects 

gue to bed friction which is related in a+ general fashion to the 
local dépth and» meen velocity, and the effect of nonlinear 


acceleration. The equation of motion is given by, 
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Figure 2.5 Definition Sketch for Exner’s €1925)° Mode. 
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Here, the KU approximates frictional effects, g, refers to. the 
acceleration due-Lo the body forces and #2 woo che pressure 
gradiant. 


The change in bed height is given by Exner's erosion 


equation, 

ony ee LP oun 

at an 3x C2 tS) 
where, an ts hemi xner's erosion scoerrilcetent. 


Bquvetion* 2. Do indicaves:*thae erosion ‘occurs it. the eilow 
velocity increases (and deposition occurs if it decreases) in 
the “downstream direction. 

Exner's: model is’ restricted to processés in which the flow 
depth is much smaller than the wavelength of the periodic bed- 
wave.” The model ¢annot’ predict the sformation of bedforms , 
but’ provides’ a-grossly correct picture “of ‘dune migrationvand 
the formation of ‘steep lee fronts: 

Gradowezyk's (1968) hydraulic model delineates wave 
propagation “in one’ dimensional -erodible ‘ped -ehannels "us ing 
shallow-water wave theory for the fluid and a continuity equa- 
tion for the bed sediment. Gradowczyk includes now-steady terms 
in the vertically integrated forms of the momentum and continuity 
equations, and indicates two 'dynamic' or surface dominated 


waves in addition to the 'kinematic' wave dominated by erosive 
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processes at thetbed:» These dynamic waves are in fact the usual 
Shallow-water waves. | 

Gradowczyk's model predicts the steepening of the down- 
stream faces of the dunes. The kinematic wave theory helps 
describe the non-linear development of bedforms, showing the 
interaction between faster-moving bedwaves and the slower- 
moving ones. Gradowczyk's hypothesis interprets the maximum 
growth conditions of bedwaves as that corresponding to the 
coincidence of crest and trough, as predicted by the kinematic 
wave theory. 

Engelund and Hansen (1966) and Callander (1969) have 
generalized the vertically integrated model to include three 
damensional flows. ‘Here the local bed stress and transport. rate 
1S assumed to, have the same direction asethe velocity vector, 
which is uniform from the bed to the: water surface’... Nonsteady 
terms are included in the basic equations, as in the work of 
Gradowezyk but their icontributions are néglected because they 
can be shown to be negligible for sediment-—induced instability. 
“The general nature of transition between dunes and antidunes 
is-understood through the model of Engelund nese at The 
detailed predictions of the mechanics of such bedforms is not 
possibile through these models since the ynydraulic model does 
not take suspended load into account. (The stability boundaries 


are critically dependent on the balance between the bed and 


suspended loads.) 


23 


2 , i _ ‘ 
ay ave o 3) ei os ae ‘ 
vay sit JORD ME. Ste eavew SLaeaINe Lekelie . B alah ? 


a al 
St aii 
tog * ; 7 : , 


R. 
~ipob @n8 tema ivecsass: on Sr Lehi 35 
ryiet eepsad S¥ee ST tRmens oat siesta ae ot Bf 


no RnIwoms *, Emini psa 7a inp tto Levsb? ABvti he «dog 

: : : ie. 
-7oeVole” ai bigs LSVEWDOY Rt a ae faye ae 
Mislxsoy e448. erengissayD efted qu Jeieeabncbans “4 
oh WribAcieseb* segs. Be 390 e068 ao en be 


why Od? wdebetobbaua ee wigs) bas PaePe t02 


oon 


ved (RG8L) pohnaide hae rao@E) noma: Sie 


i, oe 


wari ehdfoat: of) taba bes Shinada tes pase okt. are 
rte Hee Rite, eno we ‘bed De: ie, as Sah « sees 
etorlgsy viporev send eh Hiesaetetay ‘amg 48 ine 
sent aitols Sb Ritie heron pa et vad aia my is ey 
to sng Oty ph on ,eoabtan 7h 2 Le edd nt babu lieth, one ; a 
Se ee re |. bedieiged ‘ars “Bp th geen} hosid Sige SRsoWOROyD 
IRL Ces eat b Soubith =a: teahbed yh9: alt ts aot Sd 32 Cha Sd AD: 


dA 


camhidns Bve eanub coawist nofelenade tio Swder LerronoR adv. 
sal ene, gS Dillon a Tebbgmdsday! its Ore ite eodigaibms at rf) 
Sars 

tom at, Eavvcrrt Sc tsps to) edineddem sit de ooo aebioat cottased ie 


5 : m ; 

paob, fohom offiveniys si! Ponls «iebam 2nods. ceueges at ateseg > a 

i , x 4 1d 

aettabpared VILA b ae ent) ¥tivoons of 


; 
o) - 


m1 bso bshrteneys eset ton 


Dee Red si aisewted SoMB IED, att 6 titsDive noo. VLTROLI ine Sag ~ 


; | (. 2080f heabivegeag 


a7 : . | ae 

© | : ; ot 
te ri sae Vay ae 
oe ne at by ve ry 4 Sy - a . fe 4 _ Zz i Ke — o> : * fi ’ 


24 


eo SoS Potential Flow Models 

In studies based on potential flow mode iss the assumption 
of an irrotational flow of an ideal incompressible fluid sub- 
yeeu boche tnmiuence: of gravity is used. 

Although the potential flow approach was applied to explain 
bedform mechanics by Anderson (1953), it was Kennedy (1963) who 
farst used 10 dncstability analysis. 

The governing equation in Kennedy's potential flow model 


1S; 


Vor = a0 (2.16) 
where, ¢ is the potential function. 
Referring to Figure 2.6, the boundary conditions at the 


free surface are given by, 
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Equation 2.17 is the kinematic boundary condition and equation 
2.18 expresses the dynamic condition of vanishing pressure at 
therfreessurfiacé. 

The boundary conditions at the fluid-bed interface are 


given by, 


y an + an = aod On <2 = =d/ (2.19) 
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Potential Flow Model 
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Equation 2.19 expressed the kinematic condition at the bed and 
equation 2.20 is the; sediment continuity equation, where 

G(x,t) is the local rate of sediment transport per unit width 
bonwetent. and! By, as) the ‘bulky specific, weight of the sediment, in 
tHe! wb ed:. 

Kennedy's investigation concludes that the type of bedform 
and the wavelength of the bed features depend upon the Froude 
number, the depth of flow and an undetermined distance 6, by 
Wiech; the: local sediment transport, rate Lags,the, local velocity. 
A modification to the interpretation of the lag distance 6 was 
introduced by Hayashi (1970) through the addition of the 
assumption of dependance of the local transport rate on the bed 
slope as well as on the velocity. Recently Shirasuna (1973) 
propounds thay the mocion wot the erodiblee bed is composed of 
the translation of 4 sandwave;, whose travelling velocity is 
far smaller than the flow above 1G and the translation.of a 
thin surface sediment layer which flows as a fluid. Based on 
those facts Shirasuna has formed an analytical modél in which 
tne motion of theofluid and thevpiowine surtace layer atvotne 
bed are expressed as a two-layer potential flow of different 
densities. The sandwaves are treated as internal wave of small 
amplitude occurring at the interface between two potential layers. 
Although potential flow models outline the skeleton of the 
stability problem, they provide an incomplete view of the actual 
flow situations. The two severe limitations of these models 
are: 

Los There tstio, accouncyior bed shear stress. 


2. The predicted velocity variation with depth is unrealis- 
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2.6.3.4 Rotational Models 

Rotational models are real fluid models which combine ‘the 
essential features of hydraulic or vertically integrated models 
aie DOLeH Liat models, * The fimse sienifieant: contribution ain 
this respect was made by Engelund (1970), who introduced most 
of the physical processes thought to be relevant to bed stability 
into an analytical framework. 

The governing equation in Engelund's two dimensional rota- 
Crone. model js the vorticity transport equation, 


dw 
dt 


= vV?w (2..21)) 
where, w is the vorticity of a two-dimensional rate 
Both suspended and bed load are treated in the analysis. 

With the single limitation of being confined to a two- 
dimensional flow, Engelund's model retains many advantages of 
ther hydraulic models, accounting, for friction and predicting 
bed shear. It also retains an essential feature of the poten- 
tial flow model by accounting for. the role of free surface for 
short waves. When suspended load was alone considered Engelund 
found that upper range instability (antidunes) was obtained, 
while the lower range remained stable. When the bed load was 
introduced, instability occured in the lower range as well, 
taking the form ef dunes. Withvecoarse sediments, the upper 
range instability disappeared, but the antidune range was little 


ar fected with stine sediments, for. which tne bed Load is less 


important. 
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The results of the rotational model, no matter how successful 
they may be in predicting the mechanics of bedforms, apply mainly 
to sinusoidal long-crested bed waves. While these may be bed- 
forms that develop first, they are not very like the asymmetric 
randomly placed, short-crested features usually found on stream- 


beds. 
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CHAPTER 3 


STABILITY ANALYSIS 


3.1 THE THEORETICAL MODEL STREAM 


The physical processes occurring in an actual river are 


extremelyodiffaicult: tosexpress in mathematical’ terms. <A 


theoretical model river based on the following assumptions was 


eoncel vedgtoefacmlitate the- analysis of fluvial instability. 


ante 


The model river is straight and has vertical friction- 
less banks. 

The model river is free of bedforms. 

The suspended sediment concentration and all other 
associated quantities stay constant along the longi- 
tudinal direction of the river. 

The unperturbed suspended sediment concentration stays 
constant sin che tateray direction: 

Vertical eddy diffusivaty is constant, in the vertical 


direction, (Engelund 1969) and is given by, 


E sir077LUgD 


The diffusivity in the lateral direction is a constant 
multiple oft that ing the vertical directicn: 
Lateral movement of a bedload is neglected. That is, 


the equation of conservation of bed sediment does not 
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take into account that portion of the sediment load 
that moves in the lateral direction in a more or less 
close proximity to bed, i.e. the lateral bed.load. 

8. The velocity vector of suspended sediment is assumed 
to be equat to the fluid flow velocity vector plus the 
fall velocity vector. 

oO. . in-itent of assumptiom 3 and 8; if the velocity vector 
of suspended sediment is assumed constant, the erosion 
rate at the bed may be assumed “constant in the lateral 


direction. 


3.2 FORMULATION OF THE MATHEMATICAL MODEL 


3.2.1 Mass-Conservation Equation for Suspended Load 


With reference to Figure 3.1, the mass-conservation equation 


for suspended sediment for an arbitrary control volume in the 


model river can be written apwa flux balance equation in the form, 


se 7 oe f ee 1 WV eh & BME aah y Bley 5g] C3 ai 
where, c is the suspended sediment concentration by volume 
AS is the fall velocity of the sediment 
W and V are the mean Fluid flow velocities in the upward 
direction and lateral (y) direction respectively. 
Herein, by (vy) and Ey) are the eddy diffusivity coefficients in 


the lateral and vertical directions respectively. 


From assumptions 5 and 6, 


Sara Ort Us (y)D(y); 
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where B is a constant. 
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Figure 3.1 Definition Sketch for Mass-Conservation 
Equation for Suspended Load 


SE 


HoLIEvss2n09- gach wor fiosaks notsipiMad - 2.6 
. GO). Bsud Pebasqers so} Mo ttevupd 


32 


Equation 3.1 ean now be written in a modified form as: 


dc oer a. d70c d70c ac 0& 
ee eee ee ae age ly ay) | ay, | 3-8 


Here, We has been assumed to be constant and the fluid continuity 


equation for the mean flow, 


has been used to simplify the convective terms. 


Sucwc  HOrmulation of the. Boundary .Conditions 

Sacro eWavertouriace Boundary Equation 

Assuming that) in the vertical (z) direction z = 0 at the 
water surface, the equation for the water surface boundary can 
be given as, 

SOW) paee = 6 1) 
Ss OZ ise 
Z=0 

The equation indicates that the vertical flux across the water 


Ssurhace: 1S zero. 


3.c2ecre 1bed Boundary Bauations 
The bed boundary is defined as the surface where z = -D 
in the vertical direction. At the bed two conditions apply: 
1.» Bed boundary equation sor suspended sediment 
2. Bed equation of mass conservation for bed sediment 
22 seen ebedw boundary sCOnditiom tor Suspended sediment 
if 0 is the temporal rave of deposition of suspended’ sedi— 


ments on the bed from above per unit area and EA is the rate of 
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erosion from the bed per unit area due to vertical entrainment, 


then the net upward: flux at the bed is given by the "expression, 


a(vak om yor beh? = Bp D 83-5) 


But D = apn - Wel, = _p. Since the downward convection flux 
ieee 


due to fall velocity causes deposition of suspended sediment on 
Ded. 


Herice, equation. 3:5 reduces -to, 


fs = —E (3.6) 


If E is the previously defined vertical entrainment rate 


(given by E = ak J the bed boundary equation for suspended sedi- 
V 2 
Ss 


ment is given-as, 


dc Si 
Ex 5. ME ( Sale) 


For a flat bed for unperturbed: flow, lev the concentration 
at the bed be called Ch: since the rate of upward erosion must 


be equal to rate of downward deposition, 


MeGias Were (3758) 


In equation 3.8 since the flow is unperturbed, 


Cancelling vy from both the sides of the equation 3.8, 


Ch = E (3.9) 
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An appropriate evaluation for C, can be made from equation 3.7. 


Db 


3.2.2.2.2 Equation of Conservation of Bed Sediment 


Conservation of bed sediment can be applied to provide the 


retraction, 


dc = oD) 
[ -(v, - We - 7 ee Glee alo) 
moere, p = bed porosity 
meyotcaily p is equal to. .3)-~ 5.4.4 It does not affect stability 


DouUndaries and so'can be setvequal to zero in a first analysis. 
Thus, neglecting porosity, equation (8.10) can be written as: 


. s eS 
[-«v, W)c ave ie an = $2 (308707) 


Seees oecondary, Currencs zand Lateral Momentum Transfer 

Instabilities goverrieds~bDy ssediment transport have charac-— 
terispic time scales that are mucn larger, than the characteristic 
time scales of the tludd flow itseLi. ss hws it. is) valid to treaty 
the fluid flow as time-invariant during the evolution of bed 
instability. -This-assumption—is commonly -rererred—vo-as, quasi— 
steady assumption in the literature on the subject (e.g. Engelund, 
1970). The cross-sectional velocities V and W are thus et 
generated by the instabilitysitself> - The*only other mechanism 
for generating them is that which produces secondary flow. In 
an appropriately wide channel it may be surmised that these 
currents are negligible. To-expedite the analysis, the following 
assumptions are made: . 


108 Secondary thows are neglected; V = W = 0 
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Equations 3.3, 3.4 and 3.11 become, 


ac acu ac a-¢ dc d& (2a? } 
ila or el 6\ ie ' peels ge fo 35 
BRS ee = 0 (aris) 
s OZ 
z= 0 
dec oD @3. 18) 
[-Ys2 - a5] 55 
z= -D 


DL. Secondaryecurrents- provide a mechanismifor lateral 
redistribution of momentum. In the present analysis, we shall 
ignore all lateral momentum redistribution including that due 
to secondary currents, 

la. Due tovassumption ie local bedstress in the downstream 


direction tly) can be calculated from the simple formula, 
U = pgD(y)s (3575) 
where D(y) is local depth and S is the (constant) downstream 


Water surface slope. 


3.2.4 The Perturbed Equations of Sediment Transport 
The mass balance equations are considered for slight 


perturbations about steady uniform flow conditions, i.e. 


e = of (2) +e” (yszst) 
eae 
Gin aee E 
B= Dae De Cy.) 
E = ae herein, E, is the average 


entrainment rate, assumed be be constant in the lateral direction. 
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Introducing the above perturbations into the balance equations 
and reducing the equations read as follows, 


Main equation: 


2 
ay 8% Aes on : Co 14 sis) a@ct aot (3)16:) 
ce ere wae’ ee toy pt Bs 
az az OZ oY 
Water surface boundary equation: 
de 1 Z < 
O 4. dc 1 dco} = 0 
| -¥,¢, - $6 & | z=0 7 |-v.¢ ee oon eo] 
C37) 
Bed boundary equation for suspended sediment, 
2 
-é “ + Ue cea a ekorne een e ad). 0 
On aZ O 0Z az fe) d 2 
Z 
Z-—) 
(3706) 
Equation’ of conservation of bed sediment, 
de de de 
[-v =e! + [-v Sore Aig ee Ep —3 = 
So. ie =D dz 
307 % Vike | = Po + po (2 40)) 
bo OZ dz z= -D at ot 


The zeroth order terms of these relations specify the unperturbed 
concentration, which must satisfy the equation, 


2 
de dec (3.20) 
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and the boundary conditions, 


de 
ANE Et: 29 | = Or? (3.20) 
L Sing cae Ge 7, ¥=. Q 
ap StS 
dz (+205 anh = 0 (Sa ee) 


de aD 
[-v Cue € 0 | a e fe) 
we ee 2 dz Jjz= DG dat (3.23) 


The perturbed concentration e thus satisfies the equation, 


2 
Joy MMIag Ing Vy facarln See Apres tom Tos j 
Sig ee CLs Eee —e S 3 (3.24) 
P) oy az 
and the boundary conditions, 
Ex A) i tee” at = 0 (3225) 
g Z oar Ane 
dot : ace co 
az fe) 
d ane 1 lde 
Fr a fe) i ) eke fe) ul 
-—Vc +V.D —=— + &€_D es izace = ols 
Ss d fe) dz* Oo 02 az zs D. re Coo) 


Equation 3.25 is the water surface boundary equation, equation 
3.26 is the bed boundary equation for suspended sediment and 


equation, 3.2/7 is the equation of conservation of bed sediment. 


3.2.5 Evaluation of the Unperturbed Concentration 
The unperturbed concentration cy) can be evaluated from 


equation 3.20, which satisfies the equilibrium flow conditions. 
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Integrating equation 3.20 and evaluating the constants of 
integration from the two zeroth order boundary equations 3.21 


and: 3.22 results in the expression for ¢c., 
V&D 
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Herein, E has been substituted with Cb from the relation 3.9. 


3.2.6. Evaluation of Perturbed Form for & 
Referring to assumption 5, 
ER=eO/ i Us, where U, =Vt/ 0 
From assumption 10, + = pgDS 
Heplacing@el. by win thesexpression for ¢, 
& = .07)/ Vebs. b (3.29) 
In perturbed form equation 3.29 can be written as, 
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Rearranging the relation 3.30, 
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3.2.7 The Dimensionless Sediment Transport Equations 


The relations. for ¢, and a fromvequatiions..3.28 and 3.31 


fe) 
respectively are substituted into the equations Bel ee oe Seeks 
and 3.27.” The resulting equations are expressed in dimension-— 


less «f.orm.using Do as a length scale and ie as a velocity scale. 


The resulting transformations are: 
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where, the terms with tildes on top indicate them in the 
dimensionless forms. 
After dropping the tildes for sake of convenience, the 


mean equation Tor “perturbed concentration can be written as 


follows, 
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Likewise, the water surface boundary equation is, 
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and the equation of bed sediment conservation is, 
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3.2.8. The Dispersion Equation 

In hydrodynamic theory, a generally accepted procedure to 
anatyse stability is by the assumption of sinusoidal forms for 
the perturbations. The balance equations then determine whether 
or not the amplitude of the perturbations increases or decreases. 
A growth in the amplitude indicates instability and a decay in 
the amplitude indicates stabidavy:. For ‘neutral perturbations, 
the amplitude remains unchanged, corresponding to the transi- 
tion between stable and unstable conditions. 

In the present analysis, the entire lateral confinement 
of the bed is limited fo one wavelength ‘on which a4, single 
cosine wave has been imposed. The two crests of the wave are 
assumed to simulate the physical existence of banks. This is 
Piilustrated in Figures3.2. “A eréwimein thei amplitude of the 
imposed waveform, as a Tesullvotweresion ati the trough region 
and deposition at the crests, seems to correspond ‘to the physi- 
Gal process of building! bankewemiteis inf this. fashion that one 
is lked to the viewpointfof banicktormation as a stability 
process. | 

Forms for the perturbations in deptn and concentration. that 


mathematically express the above concept are, 


i t 


ec = C,(z) cos ky e” 
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D =-D,(z) cos ky e” 


Equation 3.36 introduces the idea of variation of. perturbed 


suspended sediment concentration in a sinusoidal wave pattern 
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Simulated Stream Section 


Figure 3.2 
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im the laterat=darectton. ~Baquation 3.37 Introduces the idea of 
sinusoidal bedwave in the lateral direction. 

Here k is the dimensionless fluvial instability wave number, 
related to the dimensionless wavelength 'A' by the relation 
k = = > where A= La L being the dimensional wavelength, equal 


Do 


to width of the stream. a here, is the time rate of change of 
empimicude (C, in equation 3.36 and D, in equation 3.37)«<of the 
wave. 

The balance equations can be reduced with these forms to 


yield the following. For the main equations we have, 


oo ee 
GW2) Eh ez) = Cy (alae 3e (zie = = a “Bic 
€3-59) 
The water surface boundary becomes, 
-1 
=E.0" (2) =(@yl2) 9 30°C, D e1/ey = 0 (3.39) 


The bed boundary condition for suspended sediment becomes, 


Cees 
—E,Cy(2) + =— + F Cyd = 0 (3.40) 
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and the equation for conservation of bed sediment takes the 


form, 


—E Cy (2) = Cy (Zz) + 3 CyD - = 0 (3.42) 
Z=-l 


Here, the superscript ' refers differentiation with respect to 
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Equation 3.38 is a second order ordinary differential 
equation; the solution that satisfies the two boundary conditions, 
given as equations 3.39and 3.40, can be found by standard tech- 
niques to be, 
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Ro as défined in “equations 3-43, 3-44; 3.45 and 3.46 is given 
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A solution for a can be obtained by substituting the solution 
to the main equation (3.42) into the equation of conservation 


of bed sediment (3.41)... Upon reduction, it is found that, 
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The induced instability of bankforms can be studied ier 
vermssot Pigure «3 ope (Negative tvalhvesicof satindi¢cate stapblity 
and positive values of a indicate bank-forming instability in 
the stream-bed. The figure indicates that at high wave numbers, 
greater than 0.8, the result is a stable bed. At small wave 
numbers, less than 0.8, all values of dimensionless eddy 
amit fusdivity resultilin- ansetabuiicyecivinegvrise tol the. deposition 
of suspended sediment along the sides of the channel and 
deepening the channel centre. This process is indicative of 
ines plient ibank) formation. 

The result that bank formation is initiated only at small 


wave number(s) k being given by 2mD, indicate that channels 
L 


with high width-depth ratios tend to build their banks. 
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CHAPTER 4 


SUMMARY AND CONCLUSION 


An analysis of the process of lateral erosion and deposi- 
tion of a bed of loose sediment due to flowing water has been 
presented to explain the geophysical process of bank formation 
in straight alluvial channels. 

The analysis is based upon a consideration of mass balance 
for suspended sediment. The model assumes the absence of 
secondary currents affecting the deposition-erosion process. 
The mass balance equation for suspended sediment is solved 
subject to appropriate boundary conditions at the water surface 
and bed. The solution for perturbed suspended sediment 
eoncentration is applied to -the equation, of conservation of 
bed sediment to evaluate the rate of growth of the amplitude 
of a sinusoidal bank-like disturbance imposed on the bed. A 
region of amplitude growth, or bank-forming instability, has 
been identified. Although when the amplitude becomes finite, 
the linearized stability theory no longer holds, no attempt has 
been made to include non-linear effects. 

Some results of this study are noted below: 

1; “Banks can form in. rectangular channels asa result of 

an instability phenomenon. A small bank-like distur- 
bance on an otherwise initially flat bed can initiate 


deposdeion on, the sides and scour in the centre, 
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causing initial disturbance to increase in amplitude. The 
increased amplitude of the disturbance enhances the rate 

of differential scour and deposition; presumably this process 
continues, until the banks and channel attain some stable 
Stave characteristic of some natural rivers. 

An anveresuing osopecc, Of the réesuits7or the analysis can 
explain the Tack of Dankromms in laboratory PT LUmMe Sia. cuEe 
results indicate that only at low depth-width ratios does 
bank Formation result. The typical depth-width ratios of 
fiumes beine Nigh, instaoiiicy is not expected in Laboratory 


flumes and hence no bankforms are initiated. 


—),. — (Are 
‘ae 
io.® | tae 


] 
ia 7 , rer! 
= 


on abi BE Loge ut canis om SEE Shei 
sie4 of? Eadnsdne soapdu yee abe 
stshong ebay ye same oe cep 
efdar’s sims absiia refine aOR raat 

7 Se pate t eeintion aoa” 76 5 

HES teycans "qqise te “atitibas apy Canela 

oft ~ Say m ee Af Sacer wited: eae: 
sa20 ee Te idhieninger wel te! ita odieae) 8 

Ww #6 base (RG eh wadtid qa: Sapgays a ren 
veotescdal at Ketcanae: s6erae taba, ab 
batsre bet Pir Bid an @bnett ne 

om ig are mat | 7 hes ae 7 


© 
7 } i = 4 
> i a a | Jy @ 
- bs : cs 7 ; 
: f 7 
ie om Ul) wens aig + taney 


j alk os i. iS 
pee ed i hi 


* 7 
4 * j (y al @ § ; 2 seat wQWA ) >) 
- 7 A 
1 
Uj : 
el ; ' % : 
Zi ‘ 
’ 
7 
; : : 
is f. 
( =, ral 
Pe: 7 y ; 
2 [> 
ita 
a 4 
5 & 
& 
3 
: bp : is 
; a, 
™, 
: Ve 
= : : 
7 . . fe 
o . - +, . = ! 


REFERENCES 


Anderson, A. G. (1953): The Characteristics of Sediment Waves 
formed by Flow in Open Channels, Proc. Third Midwestern 
Conférence... Fluid Mech., pp., 379~395. 


Apia. git. ate me he Rumer. Jr. (1970)<. Diffusion of ‘Sedi- 
ment in Developing Flow, Proc. off the, American soc. of 
Cy meters rr. OD. fUg= bes « 


Bagnold, R.A., (1955): Some Flume Experiments on Large Grains 
but Little Denser than the Toansportine Fluid, and ‘Their 
Implications, Proceedings Paper No. 6041, Institute of Civ. 
Baers...” wondcns pp. L6H 


Beverage, J. P. and J. K. Bulberlson, (1964): Hyperconcentration 
of Suspended Sediment, Proceedings of the American Society 
of Civil Engineers, HY6. 


Callender he A Loos): instep yy and River Channéle. 3. 
Fluid Mech. Vol. 36, Part 3, pp. 465-480. 


Binscerd.. iH. Au, WG. Anderson wand doo N. Johnson, .ClO4O is. oA 
Distinction between Bedload and Suspended load in Natural 
Streams, Transactions, American Geophysical Union, Vol. 


Einstein, H. A-s"and N.°'Chten (1954): Second Approximation to 
the Solution of the Suspended Load Theory, Univ. Calif. 
fast a thee RES. NO. 65. 


Engelund, F. A., and E. Hansen (1966): Investigations of Flow 
dy Alluviaw ooreams. heta Polytech. o¢and. ,sll—355. 


Engelund, F. A., and E. Hansen (1967): A Monograph on Sediment 
Transport in Alluvial Streams, Teknisk forlag, Copenhagen, 


Engelund, Par"@ro 7 Oy instability of Brodible Beds, J.) Rluid 
Mechs. Vole 32 part 2.00. weo—ol4. 


Engelunde ©. fe, and J. rredsoe (19 /a "* > Uhece “dimensional. 
Stability Analysas,.of Open-Channel. Flow over an Erodi ble 
Bed, Nordic Hydrol No. 2, pp. 93-108. 


Engelund, PAA. and J. Fredsoe (1974): Transition) from Dunes 


to Plane-bed in Alluvial Channels, Series Paper 4, Inst. 
Hydrodynamics and Hydraulic Eng., Tech. Univ. Denmark. 


48 


“ew 2iemns bad 26 Hoes ; 
mssVeorrth SEAT ones 


- 


ae ¥ rat meee (9 a ore fy 
Te Poe come est alt 59 2 


ee 
. 


~ 


3 Pe Dy #72 dy “tad ic es) upset a 
ist? fre Glue a; 
vi) Ta’ CaeE TARE. othe 


ness eave : LOBEL): neat et 


rod dat Agee 
YISiate novl>sedh eis $5 web ieegoe'a 


ob -STSernian oie re a) 
ie ge ae oh 


A 7ROLEL) iaoeaded: 
avivst Ab pear * ; es 
» Loy ¢ F060 Tepe pee nen 


soddemixoagga Enooag "tA MeP TF) 
Titay vin 21osq? Bao! Ba rae veoh ot 
phar ers 


wold * ese! doe thanks tage ) Neaneh = ra 27 a bn 
Bet Pee Baaee paca eee RTO , BME ST | Falta se. nt 


< 


(seam. 24 eae ay fs Ha 4 y faery aA salt 4 Hee Ape" i ae 
Nnskans $b Urey 2 Vast \, S#sea% % Ts tvul ta ar trogera@? mal 


Sigky . ,.shee sf bots sé vag Didier OTOL? 4%, - brie legend 
“| ey) ; HPS 2 Sk ; oa Pa T°SY “se Pot wii9eM 


es Lsaotensieih as'1tT if are acim at it aihet * .f | boutegnd 
Piciho tw fe Yeo wO.0S . Derierio aye / = Selb ci LBA Yate tase 
- of f ae o . a fe ‘ 9 5s bye rm pot fee 


chai tf host He age LE NY Sts “ aay AS .< bh cA a ‘Aibit dears 
Mg! occ 2a test’, aceninten Feet 


ae at 4g) bod-tnelS oy 

: Hineing= iar 2 ra je 10Vy oes 2 lngitebo tyl ; 
- - ne : 
\ e ' i nie 


4g 


Pie cher welawwB.41013)3 = Longitudinal, Dispersion. and Turbulent 
Mixing in Open-channel Flow, Annual Review of Fluid Mechanics, 
No. 8036, po. 59-77. 


Fleming, G. (1970): Particle Size of River Sediments, Proceed- 
ings of the American Society of Civil Engineers, HY2. 


Gradowezyk, M. H. (1968): Wave Propagation and Boundry Instabil- 
hey inveredyole Bed Channels, J. Fiutid Mech. 33, pp. 93-lh2. 


Graf, W. H. (1971): Hydraulics of Sediment Transport, McGraw 
Hill Book Company. 


Guy, He ©. sien. po mons. ond tw V.eRichardson, (1966): (Summary 
of Alluvial Channel Data from Flume Experiments, 1956-61, 
Geological Survey Professional Paper 462-1. 


Hino, M. (1963): Turbulent Flow with Suspended particles, 
PROC MAM @~uoOC., CiVvi. BEners.. Vol. So) sNo. HY4”. 


Hubpbedl,.oD... Won, ands D.. «Matej ka. (1959).:.. Investigations, of 

sediment Transportation, U.S. Geol. Survey, Water Supply 
Paper, 1476. 

Kalinske), Az A.s -and ¢., Pien. (1943): yExperiments, on Eddy- 
Diffusion ands Suspended, Material. Transportation, in Open 
Ghannels ..Trans.. Ais. Geoonysu. Uni0n,.NOl.o2s,. Noo 4. 


Kennedy, J. F. (1963): The Mechanics of Dunes and Antidunes in 
Frodible-Bed Channels, J. Fluid Mech., No. 16, pp. 521-544. 


Kennedy, R. G. (1895): The Prevention of Silting in Irrigation 
Canals. Min, Proce; inett 0 tet vs teiners.. Vol Cx bx. 


Lane, H. (1945): ‘"Hydrodynamics” (1st, Edn., 1879), Dover, 
New York. 


Lane, BE. W. (1941): Engineering Calculations of Suspended 
Sediment, rans. Am. Geophys... Union.) Vol. 20". 


Lene, Ho W.. (1953):9 Progress Report ion Studies on the Design 
of Stable channels of whe Bureau of Reclamation, Proc. sAm-. 
Soc. Civite engineers, Vols >. 


Lane, 1B. W. (1955):" Design of Stable Channels, Transactions, 
ASCE. Vols 1120) Paper No. 2776, pp. 1234-1260. 


Nordin. Ce Paqir.,eCl903)" sediment Transport, in Alluvial 
Channels, A Preliminary Study of Sediment Transport Para- 
meters, Rio Puerco Near Bernardo, New Mexico, Professional 
Paper 462-C, Geo. Survey, U.S. Dept. of the Interior, 
Washington, D.C. 


tastier Be! sabeaaas 
Sy Ket Pie off 


I en md ¥ 
9 Fee yeh n= To 


Ove 7. = Ler 
2 VION ot Dts. Totes 
: S oi Vt%! , 
Met | et . Sknpifames a & Dh “ 
¥ - Tae * oa 
ve J ars rE Ge f% v ankies 
“ie 


UNS (BOT) éemotite a ee bitie’ .. Beer 

‘€ pSvNlemiTeo ee es oe 

TG gegge rerutessiont hoa 
-esieti+ae bee AtLw shart. deStodayt 
(| , a is Isy™ rs rere 

fo etn het a te ac t REL) shite 
Lage “Tex eCov iG faae. mb 

“YB . tie atom sage § 7 (Rhy 


See 


Ada, ie 


ae 


neqO’ at. aeitea% : pinstiha teats 
* .of 0S: opay pi ice 5 tes 
. Slit Gre 2erng. 29 sa big make eit, 3 
: a AT a a + seh 7 eae 
i é a : ‘ i 
ain meite: te nett ya: | ‘iano it} GA horace 
Bids . gna + Ly t 154 gee oo , eters) Cy 
) ' - “3 le } 
ule’ faith ha Sr 
oe ; ofthe Sneha F imeiere aba ab sae om. al um 
; By 
hebreggPecto  anihy alieaired < Sn beg enitera bby oi ao ead: 
- 03 . gost at avy OAS ot. athe! SeeeLEeR Ey. 
Cc j 
a2 wat sé bby” ro Tas} Speeaey ot eeewy ist a anel ~ 
» ET WSS MOL! Smsi sak’ 4 Petes le 22%. 8 & to 
3 an Low Tae id Lives - 398 y i 
peaottcekerets SUITS BO tete Fo. ei MRE Me Se | sities Oe 
P feat Ee os . jo OL Sreaet OS tev gaOZA 
7 ebb em - a ts ne - 
Dc eeNGe tn NS Pog ua 2 ipumese heuer Pe ei me sibrot | , 
+2 Tt cena eres Shestne® Yo vite week Pl ord. Ah : Agee a, 
: SN Os SE6 20 “ie cos 4oM wali i e*paet gatau? of re pieten 


gror re Ceci 


“BOS Yo tae . 


‘ y Vite. a! 20 


Parker, G. N. (1976): On the Cause and Characteristic Scales 
of Meandering and Braiding in Rivers,J. Fluid Mech., Vol. 76, 
Part 3, pp: 457-480. 


Reynolds, A. J. (1965): Waves on the Erodible Bed of an Open- 
Vanier. Jo Orin nuULanMech ss VOL) 227 erart. 1, pp. 13-1337 


Reuaki Vis Ae baGloOoo): Bede Forms in) Alluvial Channels, J. of 
Plus Mechs Vols "26, Part 3, po. 507-514. 


Reynolds. Wd.) (19 (6) =) Decade's Investigation of the Stability 


of Erodible Stream Beds, Nordic Hydrology, No. 7, pp. 161-180. 


Rouse, H. (1964): Sediment Transportation Mechanics: Suspension 
of Sediment (Progress Report of Task Committee): A 
DISCUS S17 On eroCs Ame WSOC. 1 Clive Eners. s) Vol. 90, No. Yd: 


Sedimentation Engineering (1975): ASCE-Manuals and Reports on 
Engineering Practice - No. 54 


Shiracsuna. ©. (19/3): Formatdom of sand) Waves, Proc. 15th 
Congress int. “Assn. tom hydraulic Research.4., Paper ALS, 
pp. 107-114. 


Simons, De. Bs; HtOVo Richerdson. pand?) Wa... Lauens id. (1963): 
scudies: of Flow*in Alluvial ‘Channels, Some Effects of Fine 
Sediment on Flow Phenomena; Water Supply Paper 498, Geo. 
survey, U.S.) -Dept..) Oiebme Tinterior.,  Weshimeton, D.C. 


Vanoni, V. A. (1975): Sedimentation Engineering, ASCE Report 
No. 54. 


; at hoy vot a 


Gl 


“1900 26 %o bed aielatoiyl 5% 0 
CREEL 8 i eco GE. 5 


tah , etanngd®. nin ome 
| Phone ae 


adage. eis te avtesattaewuat a er 
L-TOF' gg) soe Pyeeons a eErow 


aleneiens . i fein = & ie age gerald.’ 2 teats aoe 
Roe <a gif Saiki eet to, pede Pa ae 
Ceca ee ae non 1 ego Levto ore 


To sas err ciate wie on i win 


cL +POa! . bet . eae te ro 1c 
Che ASQee + ety skeet ee Os in © 


it «tee ogsibnae pene: = 
ear rT el sda ons 


ul 


me ek 


Lat 
iu “ 
Tae 
ts 
| 
ts 
a 
& 
a i 
' 
, or 
! re 
“ f 4 
-_ b ry 
a4 x 
~ 
' , 
7) V 
‘Gee Cy eet ip ae PAP a > r , 
i a wo ‘ 7 #2 \ ' 
» AO Pa . x 4 i] ah) : : 
fee 7 we +y My by ane f ? et / } 
Bey VLSb yy WOOLY ee, et. ea) Mat Tipee bs Ain 
Als eae ./" x cs: yy "i ae oh CN 7 
A W960 uu ; “dh rE " 1 \ a ee de 
~ : ae) " 2 4 a 
; mt a i 
on \ : iy ‘i 
nt: ae 
NM Rce 
‘ wa 
“¥ Le 
Q i x ; 
aie: 
ry ye 
, ‘ ek en Fe 
; J ne 
oy 
Tey 
y ? Soy’, 
Vito 9S Bite 
i : ei; ‘ 
et che tue ret oe 
y ; i re 45 
g Wi ide ae 
i F 3 yy 
cy +1 
. , 
é Ls i - 
' 
3 7 
7 
.) 
Ye 
& 


‘ ee 


: an: mu PO we ea ee | Pree 


re 


co 


b als i 
Ree 
See A, F 7 


y Re, 
ot gee 
“we 


f 
tery 
t 

i ie 
" 


ie 


B30190 


